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Due to the space and time dependence of the wave function in the time dependent Schrodinger 
equation, different boundary conditions are possible. The equation is usually solved as an "ini- 
tial value problem", by fixing the value of the wave function in all space at a given instant. We 
compare this standard approach to "source boundary conditions" that fix the wave at all times in 
a given region, in particular at a point in one dimension. In contrast to the well-known physical 
interpretation of the initial-value-problem approach, the interpretation of the source approach has 
remained unclear, since it introduces negative energy components, even for "free motion", and a 
time-dependent norm. This work provides physical meaning to the source method by finding the 
link with equivalent initial value problems. 
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I. INTRODUCTION 



The time dependent Schrodinger equation is commonly solved as an "initial value problem" . In one dimension this 
means that the wave function ip(x, t) is fixed for all a: at a reference time, usually t = 0. However, it is also possible to 
use different boundary conditions by fixing ip(x,t) at a space point, say x = 0, at all times. This "source" approach 
has been followed to study the arrival time in quantum mechanics [[I], or characteristic propagation velocities and 
times ("tunneling times") in evanescent conditions |0-@|. In all these applications the wave function is assumed to 
vanish for x > and t < 0, 

tp s (x,t) =0, x > 0, t<0, (1) 

the potential is zero or constant for x > and t > 0, and the solution of the Schrodinger equation is sought for the 
positive half-line and positive times, ip s (x > 0, t > 0). We shall also limit the present discussion to these premises, and 
use the subscript s to denote the solutions that obey Eq. ([!]), computed for the space-time domain D + = {x > 0;t > 0}. 
They have been invariably written, up to constant factors and notational changes, as 

M*, t) = ria f" dE <? xv+,h - iEt,h Xs{E), (2) 

J — oo 

where 

1 f°° 

^ E ) = TTH dt^ s {x = Q,t)e iEt ' h (3) 

J — oo 

is the energy Fourier transform of the given "source signal" ip s (x — 0, t) and p + = (2mE) 1 ^ 2 , with the branch cut 
taken along the negative imaginary axis of E. Because of Eq. (Q), Xs(E) is an analytical function of E for ImE 1 > 0; 
however, this does not guarantee the absence of singularities of Xs(E) for real E, additional to the branch point. In 
fact, the integral over E is to be understood as going above the real axis; alternatively, E — > E + iO is to be substituted 
in the expressions above. 

It is easy to check by substitution that Eq. (S) satisfies the Schrodinger equation. Actually Allcock derived Eq. (j|) 
as the most general form compatible with Eq. ([!]) M. However he did not analyze its connection with initial value 
problems, or the meaning of the negative energies, so its physical interpretation has remained unclear. What are the 
negative energy components? How should such a wave be normalized? In time instead of position? Indeed, should it 
be normalized? Is the source to be interpreted as a multi-particle source? Except for the detailed analysis by Allcock, 
most authors have used the source boundary conditions without discussion, or as a simplifying approach whose actual 
physical content was questionable, as reflected, for example, in Moretti's words 11 The model is surely more suited 
in the electromagnetic case; for a quantum particle, it seems more satisfactory a localization in a given region, at a 
given time, rather than a "creation" during a given time interval, at a fixed position. [....] (The model) is chosen by 
reasons of simpler mathematics" . 

A "source signal" such as B H 
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ip s (x = 0,t) = e- iuot Q{t) 



(4) 



leads to an infinite norm as t — > oo for ujq > 0. An interpretation of such an infinite norm might only be possible by 
invoking a multi-particle approach which we shall not pursue here. Instead, if we keep within the one particle picture, 
it becomes clear that not all "signals" ip s (x = 0,i) may be allowed physically. This does not mean that the results 
obtained with Eq. (|4|) are meaningless. This signal may be an idealized approximation during a certain time interval 
of an actual signal, and it is also an elementary component of an arbitrary physical signal. Our basic philosophy here 
is that the physical interpretation of Eq. (|[) must rest on a process where a wave packet located in the left half-line 
is liberated at t = and crosses, in general only partially, to the right half-line. This was the point of view of Allcock 
too [0, although our scope will differ from his work. There must be a link between a physically valid ip s and some 
initial value problem, and the objective of this paper is to find and discuss such a link and the meaning of the negative 
energies in Eq. (^). The consequences on the possibility to define time-of-arrival distributions in quantum mechanics 
are also examined. 



II. EQUIVALENCE BETWEEN INITIAL- VALUE AND SOURCE BOUNDARY CONDITIONS 

A. Free motion on the full line 



We shall solve in this section the initial value problem corresponding to a wave restricted fully to the left half-line 
up to t = 0, in the interval [a, b], with a < b < 0, by an appropriate time-dependent potential that vanishes for t > 0. 
The resulting wave ipi V is only free from t = onwards. It is also useful to introduce a related wave ipf which evolves 
freely at all times, positive and negative, and is equal to ipi V for t > 0, 

ip f (x,t)=f dx'{x\U(t,t' = 0)\x')ip iv (x',t' = 0), Vt. (5) 

J a 

Here U(t,t') = ex-p[—iHo(t — t')/h] is the unitary evolution operator for the free-motion Hamiltonian Hq, 

™ lV2 



(x\U(t,t')\x'[ 



e im(x-x') 2 /2h(t-t') ^ ^ 



_ih(t-t')_ 

By setting the following source signal, 

i) s {x = Q,t)=<il Hv {x = 0,t)Q{t)=i} f {x = 0,t)Q{t), (7) 

the corresponding ip s {x, t) is obtained via Eq. (j^). We shall now establish the equality between ip s and ipf for positive 
times and positions when the "boundary condition" of Eq.(Q) is imposed on ip s . ipf(x,t) may also be written as a 
momentum integral, 



where E = p 2 /2m, and 



M*>t) ^ TT/2 / dpe^e-^/^fip), (8) 



i>M = ^ f dxe-^/ n ^ f {x,t = Q). (9) 



It can be separated into positive and negative momentum components, ipf(x, t) = ipf_ + (x, t) + ipf^(x, t), where 

V-/,±0M) = tytj f°° dpe^ h e- lEt ^ f (pM±p). (10) 

" J — OO 

The positive momentum integral is now rewritten in terms of the energy variable E — p 2 /2m, 

^ f ,+ (x,t) = W 2] q dEe^ n e- iEt / n (H) MV^E) . (11) 

To put ipf,— in a similar form, note that, because of the restriction of the initial state within the left half-line at t = 0, 
its momentum representation ipf(p), Eq. ([)]), is analytical in the complex plane p, and decays to zero as \p\ — > oo in 
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the upper half-plane. These two properties and the decaying behavior of the exponentials of Eq. (m in the second 
quadrant allow to substitute the integral along the negative real axis by an integral along the positive imaginary axis 
(there are no poles and the integral along a large arc at infinity is zero), 



1pf-(x,t) 



h 1 / 2 



dje' 



i>f(i~f) x>0,t>0, 



(12) 



where we have used p = 17 (7 > 0), and ipf(ij) is defined through Eq. (||). Changing now to the energy variable 
E = -7 2 /(2m), 



</7,-0M) = - 



h 1 / 2 



dE e ix P+/ h e ~ iEt / R 



-2E 



1/2 



i(-2mE) 



1/2 



x > 0,t > 0. 



(13) 



Adding Eq. ([11]) and Eq. ( |13[) we can now write ipf(x, t) for x > 0, t > in the form of the right hand side of Eq. (|2|), 

M*,t) = J^J dEe^^ h e- lEt / h Xf (E), x>0,t>0, (14) 



where 



Xf (E) = (m/2E) l/2 i) f \{2mE) 1/2 



(15) 



As before, the cut defining the square root branch is taken at the negative imaginary axis of E, so that the square 
roots are positive imaginary numbers for negative energies. It is also useful to write Eq. (Q) as an integral in the 
momentum plane, 



ip f (x,t) 
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dp e lpx/h e- lEt/n i}> f (p) , x > 0, t > 0, 



(16) 



where the contour is an L-shaped line that goes downwards along the positive imaginary axis and rightwards along 
the positive real axis. Let us now define the function cf> by the integral in Eq. (|l6|) without any restriction on x or t, 



<j>(x,t) 



1 



h 1 / 2 



dpe lpx/h e- lEt/h jj f (p), 



(17) 



and analyze its behavior for x > and t < 0. The contour L can be closed by a large arc in the first quadrant which 
does not contribute to the integral. By Cauchy's theorem, <p(x, t) is zero for x > and t < 0. In particular, this 
means that cj>(x = 0, t) = ip. (x — 0, t). 

We may also write Eq. (|lj]) as an integral along the real energy axis, 



<p(x,t) = 

and solve for Xf(E>) by Fourier transform, 

Xf(E) 



W-l 2 



h 1 ' 2 



dEe lxp+/h e- lEt/h Xf(E), 



dte iEt / h 4>{x = 0,t). 



(18) 



(19) 



However <f>(x = 0, t) = i/j s (x = 0,t), so that Xf(E) = Xs(E), and therefore ip s (x,t) = ipf(x,t) in the space-time 
domain D + when the boundary condition specified in Eq. (|ij) is satisfied at x — 0, and the state is initially localized 
as indicated above. 

Since in many instances the given function is the signal ip s (x — 0,t) rather than its Fourier transform, it is also 
convenient to express Eq. (^) as an integral over time in terms of the signal. By inserting Eq. (0) and changing the 
order of integration, we find the appealing expression 



ijj s (x,t)= / dt'K + (t,x;t',x' = 0)ijj s (x' = 0,t')> 



(20) 



where the kernel is given by 
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i c 00 1 r 

K+(t,x:t',0) = ~ / dEe lp+x/h e lE ( t '~ t V h = — / d PP e~ lp ^ t - t ">^ 2m ^e lpx/h . (21) 

Notice that K + (t, 0; t', 0) = <5(£' — t). Let us evaluate Eq. ( pl| ) for 2; > along the L-shaped contour in the momentum 
plane. For (t — t')<0 the contour may be closed by a large arc at 00 in the hrst quadrant, which does not contribute 
to the integral. Since no pole is enclosed, the integral is zero. This explains the upper limit t in Eq. (20). 

For (t — t') > the part of the contour in the imaginary axis may be deformed into the negative real axis. This 
allows the identification of K + with the derivative of the ordinary propagator, 

/oo 
dp-E- e V 2 (t'-t)/(2mh)^ P x/h / 22 N 
-00 rvm 

= ±jL( x \u(t,t')\o) 

mi ax 

1/2 im;r 2 

xe 2 H*-*'i (t>t',x>0). 



ih(t - t'f 

A related result was obtained by Allcock on his way to Eq. (g) 0, by considering the retarded free propagator. The 
expression of Allcock associates the source wave function to a time integral of the wave function at the source point 
times the free propagator, with a double-sided spatial derivative in between. It can be shown that it is equivalent 
to Eq. (^0|). An analogous relation is also obtained in three dimensions, with similar techniques but rather different 
objectives, by Mangin-Brinet et al ||, relating the wave function in one region of spacetime with the convolution over 
that region of the initial wave function with the propagator, plus a surface term with a double-sided gradient. 

The relation found between ips and ipf, and Eq. ( |l5| ) are the main results of this section. They provide a physical 
interpretation of the source wave function in terms of an initial value problem when there is free motion after the 
initial state is released at t = 0. Note that the knowledge of the signal ip s (x = 0,t), or of its Fourier transform 
Xs(E), allows to recover the negative momentum part of the initial state ipi v (t — 0) for the corresponding initial value 
problem only insofar as the amplitude ipf(p), known in principle along the positive real and imaginary axis, can be 
analytically continued into the negative real axis. 

In this way, the normalizability required of ijjf (or ipi V , for that matter), is translated into a property of tp s , which 
thus acquires the usual meaning of a wave-function in quantum mechanics, that is, as a probability amplitude. It 
should be observed that not all functions of the form of Eq. (||) can have this meaning: only those that are indeed 
related to a given tpi V or ijjf can be thus understood. 

We may now see the origin of the negative energies in ip s : the contribution of the evanescent waves in the source 
solution, Eq. (0), is exactly equal to the contribution of negative momenta in the corresponding freely moving wave 
packet tpf, Eq. (ra). The change of variable and the shift in the path of integration lead to what could be termed 
apparent evanescent waves, which in fact can be understood as an artifact of the Fourier-Laplace transformation in 
time. It may be surprising from a classical perspective that such a negative momentum contribution exists at positive 
times and positions, considering that the wave packet is entirely localized on the left at t = 0. In quantum mechanics, 
however, the negative momentum (equivalently, evanescent or negative energy) contribution is always present, since 
ipf t -(p) is restricted to a half line in momentum space and therefore tpf^-(x) is necessarily different from zero from 
—00 to +00 except possibly at a finite number of points. It is true however that the influence of "0/,- a t x > 
diminishes with the distance to the origin. This is illustrated in Figures 1 and 2, that show the probability densities 
\ipf(x,t)\ 2 and \ipf,±(x,t)\ 2 versus x at two different instants. The state is initially (at t = 0) the ground state of 
an infinite well between a and b, a < b < 0, and has zero average momentum, (p) =0. Its time evolution may be 
expressed in terms of know functions, see the Appendix. Note the long tails of |^/,±(^, t = 0)| 2 beyond [a, b] in Figure 
1, and the important interference pattern between ipf,+ an d "0/,- m Figure 2. 
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FIG. 1. \t/jf(x, t)\ 2 (dotted line), \tpf y +(x,t)\ 2 (solid line), and \ipf,- (x, t)\ 2 (dashed line, indistinguishable from the solid line ) 
versus x for t = 0. The initial state is the ground state of an infinite well between a — —2.01 and b = —0.01, see Eq. ( Al). 
Here and in all figures, atomic units (a.u.) are used, m = 1, % = 1. 
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FIG. 2. Same as Fig. 1, but for t = 10. 



B. Generalization for potentials that vanish outside an interval 

The above results can be generalized for an arbitrary cut-off potential profile that vanishes outside an interval [c, d] , 
located between the initial state and the source point x = 0, i.e. such that a<b<c<d<0. The time evolution of 
the wave, ipy, for x > d is given by [jTo| . |TT[l 

Mx, V = -hhj n dp e lpX/h e- lEt/h T(p)Mp) *>d, (23) 

(to be compared with Eq. (@)), where T(p) is, for p > 0, the transmission amplitude, and its analytical continuation 
elsewhere; similarly to Eq. (^), i>b{p) = ^ _1 ^ 2 / a dx e~ ipx / h ipb{x, t = 0). The contour fl goes from — oo to oo. If there 
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are bound states it passes above the corresponding poles of T(p) in the positive imaginary axis. T(p) is analytical 
in the upper half plane except in these points. Moreover it tends to one as \p\ -> oo This means that the 



same manipulations done for the free motion wave function to write it in the form of ip s can now be performed, but 
instead of ip [(2m£) 1 ' 2 ] in Eq. ( |l5|) one must write T^mE) 1 / 2 ]^ [(2mi?) 1 / 2 ] in the more general case. The same 
arguments that have enabled us to identify ipf and ip s in D + when the boundary condition of Eq. (]?]) is satisfied are 
now applicable to identify ips{x,t) and ipf,(x,t) in the same domain when 

<l> s (x = 0,t) = = Q,t)®{t) (24) 

is imposed as the signal at the source point x = 0. 

In any case free motion or motion with a localized potential barrier in the full line do not easily allow to construct 
source signals dominated by negative energy components. But these are precisely the signals of interest to study 
propagation velocities and characteristic times for evanescent conditions |^|-|| ■ The way to construct them and their 
relation to initial value problems are examined in the next subsection. 



C. "Step potentials" with different asymptotic levels 

Here we shall show the connection between source wave functions and initial value problems for the simple step 
potential, 

V(x) = e(x)V , V > 0, (25) 

as well as any other cut-off "step" potentials between c and d, such that V(x > d) — Vq and V(x < c) = 0. As before, 
we assume that the initial state is localized between a and b, with a<b<c<d<0. Following the method used in 
, it is possible to write the wave function for x > 0, after lengthy manipulations flll| , as 

i> step {x, t) = j^J dpe-^ n e^ h T l (p)i, step (p) , (26) 

where ip s te P {p) = h^ 1 / 2 dxe~ lpx / h ip s t ep (x, t = 0). This is an exact expression, where p is the momentum with 
respect to the lower level of the step; E — p 2 /2m, and q is the momentum with respect to the upper level. For real 
p, it is given by 

{ -{P 2 -Pi) 1 ' 2 P<-Pa, 
Q=\i(p 2 o-p 2 ) 1/2 ~ Po <P<Po, (27) 
{(P 2 -P 2 a) 1/2 P>Po- 

Here and in the following p equals (2mVo) 1 ^ 2 . For arbitrary values of p in the complex plane, q — (p 2 — p 2 ,) 1 ^ 2 with 
the understanding that the two branch points at p = ±po are joined by the shortest branch cut, slightly displaced 
below the real axis. In particular, for a point p = £7 (7 > 0) the corresponding q is i{^ 2 +p 2 ) ) 1 ^ 2 . T l (p) is, for 
p > Po, the transmission amplitude for left incidence, and its analytical continuation elsewhere. As in the previous 
subsection, the contour f2 goes from — 00 to 00 passing above the branch cut and possible bound-state poles of T l (p) 
on the imaginary axis. 

The integral in Eq. (|2^) may also be written in the g-plane as 

i, step (x, t)^^^ e^/ h e-' Et / h ^T l (p)^ step (p), (28) 

where now p — (q 2 + Pq) 1 ^ 2 has two branch points at q = Upo joined by a branch cut. The contour C goes from —00 
to +00 passing above the part of branch cut in the positive imaginary axis and above the poles of T l (p) (Note that 
whenever there are bound-state poles they are above the branch point at q = ipo). It is possible, as in the previous 
section, to deform the contour in the second quadrant so that the final contour has an L shape, going from ioo to 
zero passing to the right of the poles, to the right of the branch cut (p > in the interval from q — ip to q = 0), and 
from zero to 00 along the positive real axis, see Figure 3, 

ipste P (x,t) =-f^j dqe^e-^'^ip^pip). (29) 



to write Eq. (23 
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FIG. 3. Integration contours C (solid line) and L (dashed line) in the (/-plane. The crosses represent branch points joined 
by a branch cut. The circles on the imaginary axis are bound states. 



Finally, we set the zero of energy at the upper level of the potential profile using the change of variable E' = 
E — Vq = q 2 /2m, and define the wave that satisfies the Schrodinger equation for this new level in terms of the wave 
function written for the lower level, ip' step (x,t) = e %Vot ^ip a t e p{x,t), 

1 roc 

' dE'e^e-^-T^p^'ip). (30) 



~ h 1 ' 2 J-oo P 

We have thus expressed ip' step in the form of the source wave function ip s . The same arguments used to find the 
equality of ipf and ip s in D + may now be applied to identify iji s and ip s tep- In particular, 

Xs (E') = ^T l (p)jj' 8tep (p), (31) 
P 

1 /2 

where p = [2m [E + Vq)] , when 

V> s (x = 0, t) = i) step {x = 0, t)0(t). (32) 
We have performed a numerical test with the simplest step potential, Eq. (p5|). In this case, 

T\ P ) = -5-. (33) 
p + q 

The state is initially localized in the lower level of the step and sent towards the step. It is the same state used for 
Figures 1 and 2, but boosted with a non zero momentum below the step threshold. Figures 4 and 5 show the real 
and imaginary parts of Xs{E') calculated by Fast Fourier Transform of the signal %jj' step (x = 0, t)Q(t), which has been 
obtained numerically by solving the time dependent Schrodinger equation with the Split Operator Method. Note the 
clear dominance of negative energy components. The curves shown are indistinguishable from the real and imaginary 
parts of the the right hand side of Eq. (^Tj) . 



7 




-0.35 

-20 -10 10 20 

E (a.u.) 

FIG. 5. Imaginary parts of Xs(E') (circles) and (p)tp' stf , p (p) (solid line) for the step potential (Vo = 5). Initial state as 
in Figure 0. 



All the connections made in this section between the momentum amplitude tjjfip) an d the energy Fourier transform 
of the source signal at x — have required a wave packet strictly localized between a and b, in order to perform 
contour deformations. It is intuitively clear, however that relaxing this condition cannot lead to a completely different 
behavior in the relevant energy region. To test this claim, we have also calculated the real and imaginary parts of 
Xs{E') by Fast Fourier Transform of ip' step (x = 0, t)Q(t), which is obtained by solving the time dependent Schrodinger 
equation for an initial Gaussian state with negligible overlap with the positive- a; region; this is compared with the 
right hand side of Eq. ( pl| ) in Figures 6 and 7. Clearly, the right hand side of Eq. ([5l]), obtained directly from the 
initial state, reproduces the Fourier transform in the important energy region, but for very negative values of E' they 
do not coincide. 
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FIG. 6. Real parts of Xs(E') (circles) and (p)ip' s te P (p) (solid line) for the step potential (Vb = 1). The initial state is a 
minimum uncertainty- product Gaussian with (a;) = —3, (p) = 1, and A x — 0.5. 
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FIG. 7. Imaginary parts of Xs(E') (circles) and (p)tp' stf , p (p) (solid line) for the step potential (Vo = 1). Initial state as 
in Figure H. 



III. APPLICATION TO THE ARRIVAL TIME 



Allcock obtained Eq. (Q) by assuming that Eq. (|l]) was satisfied, and that V(x) — for x > 0, without any further 
specification of the potential after the initial state is released jl| . But each experimental setup will involve a particular 
potential profile, and we have seen that the interpretation of the source solutions in terms of an initial value problem 
depends on the particular potential that determines the dynamics. We shall in particular choose the simplest case of 
free motion to analyze Allcock's objection to the possibility of defining an ideal (apparatus independent) arrival-time 
concept in quantum mechanics based on Eq. (]|) . This section is an elaboration of the arguments given in ju| . 

Allcock's motivation for studying quantum mechanics with sources was the fact that the expression for the wave 
function, Eq. (^), involves an integral over energy from — oo to oo. This suggested a possible way out to Pauli's 







argument against the existence of self-adjoint time operators conjugate to a semi-bounded Hamiltonian fll7| ]. 

The core of Allcock's objection rests on the identification of the total arrival probability at a point X > 0, P(oo), 
with the total amount of norm to the right of X as t — > oo, 

y* oo />oo 

P(oo)=lim/ dx\^ s {x : t)\ 2 = dE\{E\^ s )\ 2 . (34) 

t-H^Jx Jo 

Here 

(£|V^X S (£)(2£/™) 1/4 , (35) 

where, as before, the branch cut for the fourth root is set at the negative imaginary axis. For free motion the second 
equality in Eq. ( p4| ) is to be expected since (E\ip s )(2E/m) 1 ^' i , with E > 0, is nothing but the amplitude for positive 
momentum p, so that we recover from Eq. (B4I) the known result 



lim / dx\il> a (x,t)\ 2 = dp\iP f (p)\ 2 . (36) 
t^^Jx Jo 

Using two other independent routes Allcock arrives at integral expressions for a hypothetical total arrival time proba- 
bility in contradiction with Eq. (|34| ) because they have contributions from E < 0; the simplest one is dE\ (E \ip s ) \ 2 . 
He therefore concludes 11 unequivocally that an ideal concept of arrival time cannot be established for the problem with 
sources (— oo < E < oo)". In the light of the identification between ip s (x,t) and ipf(x,t) in D + made in section II, 
the flaw in Allcock's argument is that he is overlooking the contribution of negative momenta. According to classical 
prejudice one would not expect that negative momenta play any role for t > at X > 0, if the initial state is confined 
within the left half-line. But quantum mechanics works otherwise, as Figures 1 and 2 clearly show. Since the negative 
energy part of tp s (x, t) is equal to ipf_-(x, i), we may equivalently think of the contribution of the negative momenta 
to a time-of-arrival probability in terms of evanescent waves. The negative energy components penetrate the posi- 
tive coordinate region for some time before finally abandoning it. But Allcock's definition of P(oo) disregards this 
contribution entirely. Note that our criticism to the validity of P(oo) is rather independent on the theory chosen for 
defining an arrival time distribution. However, we cannot fail to point out that the e quality between tJ) s and ipf in the 
space-time domain D + allows immediately to associate Kijowski's distribution Hk |l8|,|l6| t o the source problem for 
free motion. This arrival time distribution has been extensively discussed elsewhere [|l9|-|2l|| . Similarly, for a general 
potential profile the generalization of Hk proposed in |]20|| is applicable. 



IV. DISCUSSION 



The physical meaning of source boundary conditions used for the one dimensional, time dependent, single particle 
Schrodinger equation has been clarified by showing the equivalence of the corresponding solution ip s , see Eq. (g) , with 
standard initial value problems for positive times and positions. We have first shown the connection between Eq. (|^) 
and the ordinary initial- value-problem integral expression for a free motion wave packet, and then have extended 
this result for other potential profiles, in particular for the step barriers, which allow to inject waves dominated by 
evanescent energies into the right half-line. 

Each of the potential profiles may correspond to a different experimental arrangement. Depending on the potential 
the origin of the evanescent waves in the source solution may be traced back to either negative momentum components 
(which are present in all cases, and become the only source of negative energies for free- motion), or "true" evanescent 
components, namely energy components which are already negative in the associated initial value problem, relative 
to the region of positive x. This occurs for example in the step potentials, where positive momenta with respect to 
the lower level become imaginary momenta in the upper level region. In addition to scattering contributions to the 
negative energies, we have also allowed for the possibility of bound state contributions, that may become important 
if the initial state overlaps with the normalized eigenstates of the Hamiltonian. 

Apart from elucidating the physical content of previous works that have made use of source boundary conditions, 
our results have also served to analyze critically Allcock's negative statements about the possibility to define an ideal 
probability distribution of time of arrival for quantum waves with sources. It is argued that Allcock's postulated 
expression for a total arrival probability is flawed since it ignores the quantum contribution of negative momenta at 
positive positions and times. 
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APPENDIX A: SOLVABLE EXAMPLE 



We shall work out the time dependence of a state which is at t — the ground state of an infinite well of width D 
between a and b = a + D, and evolves freely from that time on. This was the only example of source mentioned by 
Allcock J§: 



2X l/2 



ipiv(x,t = 0) = I — 1 sin[(ar- a)k w ]H(a,b), (Al) 
where k w — tt/D, and 

H(a,b) = { 1 \< X<b . (A2) 

1 otherwise 



Eq. ( Al) may also be written as a combination of four plane waves truncated at a and 6, two of them with momentum 
Pw = k w h, and the other two with momentum — p w . The time evolution of a truncated plane wave was first solved 
by Moshinsky |2^] . For an alternative derivation see [nj . The resulting wave is given by 

1 / o \ 1/^ 

where 

±p w — m(x — b)/t 

1 ' 
±p w — m(x — a)/t 



Ub± 



u a ± 



f 

/ = (!-*) f 



* \ 1/2 

mn x 



t 

w(u) = e~ u ~ eric(-iu). (A4) 



For properties of the w- function see [|23|. The momentum representation of Eq. ( Al) is also of interest to illustrate the 
properties mentioned in section II (analyticity and decay to zero in the upper half plane when the state is confined 
between two points in coordinate space), 

Oh} 1 / 2 « , p -i(p+aPu,)b/h _ -i{p+ap w )a/h 

^•(P) - -TJym E a ~ T 1 e ■ (A5) 

^U 1 ^ p + ap w 
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